Stationary entanglement in strongly coupled qubits 
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The dynamics of two superconducting flux qubits coupled to each other and to a common bath is discussed. 
We focus on the case in which the qubit-qubit coupling strength dominates over the respective qubit transition 
frequencies. We derive the master equation including collective effect by modeling the bath as ID open space 
in this ultra- strong coupling regime, and find that the coupling greatly modifies both the coherent and the in- 
coherent dynamics of the system, giving rise to qualitatively diff'erent properties. By analyzing the steady-state 
and the dynamics governed by the master equation, we show that ground state entanglement and maximum 
coherence between the two qubits can be induced by the environment alone. By employing in addition a single 
external driving field, both the entangled anti- symmetric and symmetric collective states can be populated and 
preserved with high fidelity. Similarly, entangled states can be prepared using adiabatic passage techniques us- 
ing two external fields. Our results could find applications in entangling quantum gates and quantum memories 
free from the decoherence. 

PACS numbers: 85.25.Cp, 03.67.Bg, 03.67.Pp, 42.50.Hz 



I. INTRODUCTION 

Entanglement is the primary resource for quantum compu- 
tation and information processing \ and correspondingly, effi- 
cient entanglement generation and preservation is a subject of 
intense research in all physical systems considered as possible 
realization for quantum protocols. Here, we focus on super- 
conducting qubits^'^, which in principle have the particular ad- 
vantage of great freedom in designing the system properties. 
They promise an attractive route towards a scalable quantum 
computer, such that creation and protection of entanglement in 
superconducting qubits is of particular interest. Consequently, 
entanglement generation between superconducting qubits has 
been suggested and demonstrated recently by a number of 
groups^^^. Usually, the creation of entanglement requires 
precise control of all fields interacting with the qubits. How- 
ever, in typical experiments, the inevitable decoherence and 
the imperfection of actual driving fields reduce the achievable 
fidelity for state preparation and quantum gates. 

The decoherence, which arises due to the coupling of a sys- 
tem with its environment, usually tends to destroy the quan- 
tum effects and entanglement at the heart of quantum proto- 
cols which provide the advantages over corresponding clas- 
sic algorithms. Superconducting qubits are particularly lim- 
ited in their coherence time due to the many uncontrollable 
degrees of freedom in their solid-state environment. Sev- 
eral proposals have been made to eliminate noise-induced 
decoherenceA^ii2r— . In particular, by careful engineering of 
the local environment coupling to a system, entangled pair 
states have been prepared using a dissipative processi^^i^. 
However it is difficult in experiments to isolate well the lo- 
cal environments coupled to the individual qubits from each 
other. Moreover, for a given setup, the steady state which can 
be prepared is unique. Alternatively, entanglement also has 
been created between two qubits interacting with a common 
heatbath^^'i^ 

Next to the creation, also the preservation of entangled 



states is a long-standing challenge. As a possible solution, 
it has been pointed out that ground state entanglement is in- 
trinsically immune to the decoherence. It has been discussed 
in complex system22i2^, but without addressing the inevitable 
decoherence. An alternative is the well studied steady state 
entanglement"^'^"^. You et al.^^ discussed a highly entangled 
ground state based on analyzing the coherent Hamiltonian 
without taking into account the dissipation. Certain kinds of 
ground state entanglement in superconducting circuits have 
been discussed in topological computation^^ recently. The 
possibility of ground state entanglement has also been dis- 
cussed in the situation of a vanishing separation between two 
superconducting qubit a^^i^^ . But in these references, neither 
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FIG. 1 . (Color online) The considered model system. Two supercon- 
ducting flux qubits are coupled through their mutual eff'ective induc- 
tance M and coupled to a common reservoir modeled as an LC circuit 
via the mutual inductance Mb. We consider the case of strong cou- 
pling between the qubits, in which the always-on coupling strength 
between the two qubits exceeds their respective transition frequen- 
cies. 



the coherence nor the time evolution have been considered. 
Experimentally, ground state entanglement with concurrence 
C = 0.33 between two inductively coupled flux qubits cou- 
pled to a common heat bath has been observed^. In this work, 
the always-on coupling (AOC) strength is close to the small 
transition frequency of qubits. This suggests that strong qubit- 
qubit interactions could be an interesting parameter regime for 
the creation and preservation of entanglement. This is further 
supported by related work on atoms subject to intense exter- 
nal driving fields, which can lead to a substantially modified 
incoherent dynamics of the system such as the population of 
excited bare states^^ . 

Therefore, in this work, we analyze a system of two inter- 
acting qubits coupled to a common bath. In contrast to pre- 
vious results, we focus on the case of ultrastrong coupling 
between the two qubits, with coupling strength exceeding the 
respective qubit transition frequencies. We derive the corre- 
sponding master equation, interpret it in diff'erent collective 
state bases, and analyze the steady state properties and the 
temporal evolution without external driving fields, and with 
one or two external fields. The aim of the analysis is the gen- 
eration of stationary entanglement between the two qubits. 

We find that both the incoherent and the coherent dynam- 
ics is qualitatively modified by the strong always-on-coupling. 
As a consequence, ground state entanglement can occur even 
in the absence of any driving field, with the entangled state 
populated by spontaneous emission only. By dynamically 
switching the strong qubit-qubit coupling on and off, we found 
that starting from an initially unentangled ground state of the 
two qubits, entanglement with high fidelity can be generated 
on time scales on the order of the inverse qubit transition fre- 
quency. With the help of a single external driving field, both 
symmetric and antisymmetric entangled qubit states can be 
prepared and maintained with high fidelity. Finally, using two 
external fields, entangled states can be trapped based on so- 
called dark-states or prepared via adiabatic passage methods. 
Possible applications of our scheme include nonlocal quantum 
gates and quantum memories. 



II. MASTER EQUATION 
A. Model 

Our system consists of two flux qubits coupled to each other 
through their mutual efl'ective inductance M and to a reser- 
voir of harmonic oscillators modeled as LC oscillators, as 
shown in Fig.[T] Both qubits are operated as two-level quan- 
tum systems. A crucial difl'erence to previous work^'^^ is that 
in the present work, ultra- strong coupling due to a large mu- 
tual inductance M is considered. This coupling can be engi- 
neered22i2^, and is dynamically controllable^—. In particular, 
it can be comparable tc^^'^^ or even much Iargeri22i2i^ than the 
qubits' transition frequencies. 

Each qubit is modeled as a loop interrupted by three Joseph- 
son junctions: two identical junctions characterized by the ra- 
tio of Josephson energy to charging energy Ej /E^J, and a 
third one with area smaller by a factor ^/ than the other two. 
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FIG. 2. (Color online) The transition dipole matrix element as a func- 
tion of bias ft = OJOo for ^ = 0.8, Ej/Ec = 35. 



Here, / e {1,2} labels the two qubits. Thus, the Josephson 
energies and capacitances in the /th qubit loop are given by 
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The gauge-invariant phase drops across the three junctions in 
the /th qubit are 0^^\ 0^^ and ^3^ The qubits are modeled as 
two-level system, and in pseudo-spin operator language, the 
single-qubit Hamiltonian can be written as 



H\ 



-\W 



® + 6,<rfl. 



(2) 



This model is valid for both three- or four- junction flux 
qubits^^. The transition frequencies 6/ = llf((b^^ - Oo/2) 
are tunable via individual static bias magnetic fluxes O^^ char- 
acterized by the bias ff = O^^V^o, where Oq = h/(2e) is the 
superconducting flux quantum. The persistent current through 
qubit lis Ip\ and the tunnel coupling between the two poten- 
tial wells is denoted as ti. Additionally, each qubit is driven 
by a TDMF 



^f(t) = Ai cos(ajft) , 



(3) 



with amplitude A/ and frequency ojc . The coupling between 
the two qubits is modeled using a generic Ising-type coupling 
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)cr. 
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We operate in the small vicinity of the optimum point /^^ 
0.5, and use a rotated coordinate system such that cr^ and cTx 
are efl'ectively exchanged^^'^^. The coupled qubit pair driven 
by the TDMFs can then be described by the following Hamil- 



Ionian in a bare- state basis2ii2 
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- n Y, [o-^^e-''^^' + H.c.) (kie-"^''' + c.c) 
-nj]{cTfcT'-^+H.c.)(^y^'' + c.c) 
-/.|]KVf+H.c.)(Qj>-?' + c.c) 

+ 7z(MV5 + MVf + /f.c.). 

The transition frequency uXq of qubit / is determined by 



/^? + f'- 



(5) 



(6) 



The Pauli spin matrices cr+ _ ^ of the /th qubit with ground 
state 1^/) and excited state |e/) are defined as 



af = \ei}{ei\ - \gi}{g,\ 
<r^^ = \ei}{gi\, 
<r? = |g/)<e/|. 



Introducing 
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QfJ = Ai/Siiei, g„|/<r" cos(2^® + 2;r/f)|g,, e^}/2 , (9b) 
nfj = AiMei, e^\lf^ cos(2.^« + 2;:/®)!^,, g^)l2 , (9c) 

and can be controlled by the applied TDMFs. Here, the pa- 
rameters are given by^ 



4> = (cf,f+cf,f)/2, 
0® = (<-^<'>)/2, 

the coupling strengths can be written as^ 

k=Ai{ei\(lf + iriPpj)\gi}/(2n), 
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the effective momenta are 
Ppy 
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and the dipole moments are defined as^ 
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= <e,|sin(2<^^" + 2;r/;'0lg/). 



The always-on coupling (AOC) energy arising from the mu- 
tual inductance induced qubit-qubit coupling is given byM 



J^MI^Hf. 



(13) 



This AOC strength is real in the vicinity of the optimum point 
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f^^ = 0.5. The ultra- strong coupling regime is characterized 
by / exceeding the transition frequencies co^ of the qubits. 

Interestingly, the transition dipole matrix elements can 
change their signs as the bias fluxes sweep across the opti- 
mum points, see Fig.O We define the parameter 



\di\\d2\ 



(14) 



characterizing their relative orientation with p = 1 denoting 
parallel and p = -1 anti-parallel dipole moments. 



B. Diagonalizing the Hamiltonian 



As we focus on the strong coupling regime, a perturbation 



treatmen 
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of the coupling / cannot be applied. There- 
fore, we transfer the system to a suitable dressed state basis, 
thus taking into account the AOC to all orders. The diago- 
nalization of the Hamiltonian in Eqs. ([5]) without TDMF but 
including the AOC leads to the eigenenergies Ej and corre- 
sponding eigenstates |j) (j e {1,2, 3, 4}), which can be inter- 
preted as a single four-level system with states and energies. 
Note that the states and energies are diff'erent from the usual 
collective Dicke states commonly introduced for two weakly 
interacting two-level systems^^. Instead, due to the strong 
coupling, the eigenstates and energies are given by 



with 



\l} = a\eue2} + b\gug2}, Ei=nw, 
\'2}=/3\eug2} + a\gue2), E2 = nw, 
\3} = a\eug2)-/3\gue2), E^ = -Hw , 
\4-} = -b\eue2} + a\gug2}, E4 = -Hw , 
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D = cl)o + w, m = A + w, 
D . J 
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In the following investigation, we assume two identical qubits, 
i.e., A = 0. This reduces the above expressions to w = |/|, 
a = -^, 6 = -^—, and a^ -\- b^ = 1. This assumption is 

V2 ^ V2i/r ^ 

reasonable as long as A «: |/|. In terms of the new dressed- 
state basis, the coherent Hamiltonian of our system including 



the TDMFs becomes H = Hq -\- Hj with 
^0 = / EjRj 






(17a) 



Hi =nQi {(aa + bfiXRu + 7^34) 

Hcij3 + ba)pe-''^''(Ru - R34)] e''"^'' 

-^n^2{(afi-ba)(R24-Ri3) 

Haa - bj3)pe-''^''(R24 + ^13)} ^"'"'^^' 

+ n^3(a^ - b^) (1 + e-''^'')e-''^^'Ru + H.c. , (17b) 

where the PauH operators are defined as 

^.7 = |/)01 (18) 

for /, J G {1,2,3,4}. The interaction Hamiltonian Hi arises 
from the coherent driving via the TDMFs. For simpHcity, we 
have assumed that the field frequencies and Rabi frequencies 
satisfy cOc = (^c = ^l, ^1 = ^2 = l^il = Ifel and ^3 = 
|Q^^2 1 = 1^21 1- ^^^ relative phase of the driving fields between 
the two qubits is 

021 = hrii = 27T(r2i/Ao)(ojL/oJo) . (19) 

C. Coupling to the bath 

We model the interaction of our qubit system with the envi- 
ronment via a reservoir of harmonic oscillators. The system- 
bath coupling can be described using the Jaynes-Cummings 
Hamiltonian^^*^ 

2 

k /=1 k 

+ ^ Z Z ^* (^*"«* + ^•^•) (^®^?" + H.C.) , (20) 



k /^m=l 



where r/ is the position of /th qubit. The bath oscillators have 
frequencies cok, and (alak) = rik is the average photon number 
of the ^th field mode. The thermal average number of photons 
for an oscillator of temperature T is given by 



Nth(oj) = [cxpihoj/kBT) - 1]- 



(21) 



and is assumed negligible at the frequencies relevant to our 
system, t]). and Xr are determined by the vacuum field 
and are proportional tc^ ^bA ^^<^r/2LB, where M^ is the 
qubit-bath mutual inductance in the corresponding qubit. For 
simplicity, we assume M^ = M^ = Mb. In supercon- 
ducting circuits, the qubits can be driven or biased by mi- 
crowave circuits^i22i22i^ or ID transmission lines^^"^^. Thus 
the bath can be modeled as one-dimensional open spaced, 
and an Ohmic environment with a large cutoff' ojcut is a re- 
alistic assumption for this bath^Sii. 

We assume that the coupling to the bath is weak, i.e., we 
operate in a Hamiltonian-dominated regime. Both qubits are 



operated near their optimal point, so that they are decoupled 
from low-frequency noise to a good degree, which otherwise 
would give rise to the pure dephasing. Then, applying the 
Born-Markov approximation to eliminate the batbi^s^Si^, we 
find that the bath introduces both incoherent and coherent con- 
tributions to the dynamics of system via their dipole-dipole 
interaction (DDI), which we discuss next. 



1. Bath-induced incoherent contributions 

The incoherent part, i.e., the bath-induced dissipation, can 
be described by the Liouville operators ^spP and ^tpP 
operating on the system density matrix p. The fist part 
^spp results from single-photon process proportional to rik 
in Eq. ([2Ql) . It is composed of two contributions as 



where 



^,PP = ^iP + ^2P, (22) 



^ip =^{r| 1)^12), 7?2l}p + ^{r|2>^|4>,^42}p , 

+ ^{r™,^i2,^24}p, (23a) 



^2P =^[y\l)^\3),R3l}p + ^{r|3>^|4>,^43}p 



+ ^{r™,7?i3,^34}p, 



^{k, A]p =^ (2ApA'^ - A'^Ap - pA'^a) 
^{K,AuA2}p =k(aIpAi + A]pA2) . 



(23b) 
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The rates are given in Appendix (A] The decay rate of an iso- 
late qubit is given by 

yo(oj) = sm\6)(MDlpfoj/(2nZ) . (24) 

This rate describes the damping of a qubit coupled to an ID 
open transmission line with impedance Z = ^|lr|cr by a mu- 
tual inductance Mb"^^. Here, Ir and Cr are the inductance 
and capacitance per unit length, respectively. In Eq. (l24l) , 
sin(^) = tjoj with tunneling t and the transition frequency co. 
In our case, sin(^) ^ 1 . Note that an isolated qubit decays at 
a rate 270 in our case. These results are in agreement with 
Ref. H The dipole-dipole crossing decay rate 712 evaluates 
to 



yn^co) = 7o(^) cos(Kr2i) , 



(25) 



where k = cdIv is the wave number corresponding to the tran- 
sition frequency and the wave phase velocity v = 1 / ^]lrCr. 
The eff'ective separation between the two qubits is r2i. This 
formula was already presented in our previous work"^^, and a 
similar superradiant decay rate was recently found in a donor- 
acceptor system mediated by plasmonic nano waveguide s^2i^. 
The second part of the bath-induced decoherence ^tpP is 
due to two-photon processes and can be written as 



^TPP = ^{yTP,Ru}p, 
with the relaxation rate 

yjp = 4(a^ - b^fifo -\- fu) , 



(26) 



(27) 



where the two-photon decay rate 70 is defined in AppendixlAl 
The rate fu is given by 



712 =focos(Kr2i). 



(28) 



All relaxation rates are derived in Born-Markov approxima- 
tion neglecting small Lamb shifts^?^, and are summarized in 
AppendixlAl 



2. Bath-induced coherent contributions 

The coherent part Hcddi consists of dipole-dipole energy 
shifts 



Hoos=nY,^f^''\J)U\ 



J 



(29) 



in the transition frequencies and of the bath-induced excita- 
tions Hbie- These two contributions are discussed in detail in 
Appendix [Cl 



D. Master equation 

After including the incoherent and coherent contributions, 
the dynamics of our system can be described by a master equa- 
tion as 

dp . v^ . „ ^ i 



op Y^ I 

_ = - / 2^ ojj[Rjj,p] -TVtii + Hbie.p] 



dt 



-h ^spp + ^TPP + 



(30) 



Here, we redefined the transition frequency to also include the 
dipole-dipole shifts as 



fiojj - Ej -h ficxr. 



(DDS) 



(31) 



The final term ^^ phenomenologically introduces pure de- 
phasing with rate F^ for the two qubits to the master equation 
Eq. [30l and is given by 

^^P = Yj^ {r0/2, \ej}{ej\ - \gjXgj\]p . (32) 

j 

Note that this treatment of pure dephasing is only valid for 
operation near the optimum point, where the qubit in leading 
order decouples from fluctuations in the external fields, and 
thus becomes insensitive to dephasing induced, e.g., by 1// 
noise. 



E. Model parameters 

As already mentioned, we assume identical qubits and 
equivalent TDMF parameters for the two qubits. Further- 
more, for our calculations, we assume a large two-photon de- 
cay rate 70 = O.O270 and r2i between /lo/lOOO and /lo/50 cor- 
responding to an eff'ective separation of around 10 - 200 jam 



between the two qubits22iS. Note that the separation r2i/A 
scaled by the wavelength even can be dynamically changed 
in superconducting qubits because the transition frequency oj 
is tunable via the bias /^. It also can be influenced by aux- 
iliary qubits^ ^'^^. The qubit transition frequencies are taken 
as ^0 = IOOO70, and the bath cutofl' as ojcut = 200a;o. We 
adjust the bias fluxes to operate the qubits near their optimal 
points fy = 0.5 and assume identical, but uncorrected pure 
dephasing for the two qubits. 



R Experimental realization 

We now discuss aspects of the experimental realization of 
our framework based on existing experiments and proposals. 
As our analysis focuses on the case of strong coupling be- 
tween the qubits, it is important to note that such a strong AOC 
between two flux qubits comparable or even larger than their 
transition frequencies has been realized22i22i2ii^2i^. This cou- 
pling can also be dynamically controlled in experiments^—. 
These advances open the door to the study of an entirely new 
parameter regime and motivate our work. 

Regarding the qubits themselves, in practice, two fabricated 
qubits will always be slightly difl'erent, e.g., in the transition 
frequency at the optimal point. But this efl'ect of nonidenti- 
cal qubits is to a good degree negligible since the coupling of 
interest in our work is at least comparable to the transition fre- 
quency, and thus strongly dominates over possible detunings 
of the two qubits. 

Next, we discuss the bath, and focus in particular on its ID 
nature, and the fact that we treat the qubits as coupled to a 
common bath. In our system, the superconducting qubits are 
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or 



driven by microwave fields led by microwave lines 
transmission lines^^i^iS. In both cases, the spontaneous emis- 
sion is the dominant source of energy relaxation^^^SiS^iS^, 
which results mainly from high-frequency noise. In a recent 
experiment^^, a flux qubit is coupled to a ID transmission line. 
The qubit is well isolated from the other degrees of freedom 
in the surrounding solid-state environment and behaves as a 
nearly isolated artificial atom in open space, coupled only to 
the electromagnetic fields in the ID space with long correla- 
tion length. In this work, the authors presented an analytical 
expression for the relaxation rate which is in good agreement 
with the experimental observation. Storcz and Wilhelm^ pro- 
vide two flux qubit coupled by flux transformers^ as an exam- 
ple for a situation described by a common bath. Note that both 
microwave lines and transmission lines induce an Ohmic envi- 
ronment and can be modeled as a global LC oscillatorii^i^iS^iS 
as in our analysis. Thus the above works indicate possibilities 
to implement our framework. 

We have furthermore assumed a bath of negligible tem- 
perature at the qubit transition frequencies. A bath of finite 
temperature would gives rise to an incoherent pumping to the 
excited state of system. Thus, the system would evolve into 
a mixed state in thermal equilibrium even if no driving field 
is applied. This can lead to a reduction of the entanglement 
found in some of the cases discussed below. Thermal pump- 
ing can best be avoided by using qubits with large transition 



frequency embedded in a low-temperature refrigerator. 

Next, we turn to the effects of the dipole-dipole interaction, 
which arise from the coupHng to a common bath. Due to this 
coupHng, two qubits decaying to a common bath may exhibit 
collective effect such sub- and superradiant states. Recent ex- 
periments^2i^ have observed such superradiant and subradi- 
ant states of two distant superconducting qubits coupled via 
cavity modes. The dipole-dipole interaction of the qubits in 
these two works is a result of exchange of virtual photons via 
the cavity. This interaction is also the origin of the collective 
effects. This cavity can be replaced by flux transformers^ 
or a common readout or driving circuit^'^^. Two adjacent flux 
qubits coupled by magnetic dipole-dipole interactions are also 
discussed in a recent review by Clarke and Wilhelm^ . Finally, 
the DDI between two flux qubits can also be induced by flux 
transformers^^ or a quantum bu*^^'^^ . 

Finally, we discuss decoherence. The decoherence consists 
of two contributions. The first contribution arises from the 
relaxation process, and the second is pure dephasing arising 
typically from low-frequency noise such as 1 //-noise. Our 
system is operated near the optimum point /, = 0.5 of the 
flux qubits, in order to reduce the effect of the pure dephasing. 
Then, the qubits to first order are decoupled from the dephas- 
ing fluctuations^'^"^'^^. This allows us to phenomenologically 
treat the pure dephasing via a Markovian master equation. 



III. RESULTS 
A. Discussion of master equation 

In this part we discuss some of the properties of the master 
equation (l3Ql) , and compare it to corresponding master equa- 
tions in atomic or ionic systems. We find that the implemen- 
tation in strongly coupled qubits offers a number of distinct 
features. This mainly arises from the fact that in our model, 
the superconducting qubits are not only coupled via the bath, 
but they can also be directly coupled via their mutual induc- 
tance. This coupling includes a two-photon channel, which 
is not present in typical atomic/ionic systems. Furthermore, 
we consider the ultra-strong coupling regime, which appears 
unfeasible in atomic systems 

We first note that the dissipative dynamics is strongly mod- 
ified by the strong AOC via the mutual inductances. In atomic 
systems, the relaxation is determined by the collective dynam- 
ics depending on the separation of the two particles and their 
mutual dipole orientation^. In our system, the effect of dis- 
sipation can most easily be described in the dressed state ba- 
sis shown in Fig. [3j and defined in Eq. ([T5]) . According to 
the master equation, as expected, the upper states decay to 
the lower states, such that all population resides in the ground 
state |4) in steady state in the absence of a driving field. How- 
ever due to the strong AOC, this ground state can have a rather 
complex structure. To establish this result, we transfer our 
system into the collective Dicke state basis with upper state 
\E) = 1^1, ^2) and lower state |G) = l^i, ^2) as shown in Fig.lH 
This basis is frequently used to analyze weakly coupled two- 
qubit systems. In this representation, the coupling to the en- 
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FIG. 3. (Color online) The relevant relaxation channels in the col- 
lective dressed state basis. The transition |1> -^ |4> is a two-photon 
channel. 
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FIG. 4. (Color online) The system in the collective Dicke state ba- 
sis. The black downwards arrows denote dissipation, whereas the red 
upwards arrows indicate pumping which is introduced in the collec- 
tive state basis due to relaxation in the dressed-state basis shown in 
Fig.El 



vironment not only leads to relaxation from states with higher 
energy to states with lower energy, but also effectively to 
pumping from lower to higher states, and even to some coher- 
ent interaction (see Appendix|B]). Interestingly, both damping 
and pumping emerge in every transition channel, even though 
the pumping is always smaller than the damping (\b\ < \a\). 
Thus we can conclude that already the dissipative dynamics 
can lead to the formation of states of interest, and we will 
show below that this also includes highly entangled states. 

We now turn to a discussion of the collective state decay 
rates indicated in Fig. [3l In the weak-coupling regime, these 
cooperative relaxation rates are essentially determined by the 
effective separation r2i between the two qubits and their tran- 
sition frequencies. But in our model system, the relaxation 
rates in addition crucially depend on the strong AOC between 
the qubits. As an example, we show the dependence of the 
collective decay rates on the AOC strength / in Fig. \5\ For 
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FIG. 5. (Color online) The collective decay rates yi^2 [(i) thin blue 
line], 72^4 [(ii) red line], 71^3 [(iii) green line] and 73^4 [(iv) brown 
line] as a function of the AOC coupling /. The distance is chosen as 
rn = /I0/5O. In (a), the two transition dipole moments are parallel 
(p = 1). The result for anti-parallel dipole moments (p = -1) is 
obtained by mirroring the figure at the / = axis. 




FIG. 6. (Color online) The collective decay rates 71^2 [(i) thin blue 
line], 72^4 [(ii) red line], 71^3 [(iii) green line] and 73^4 [(iv) brown 
line] as a function of the coupling / for r = /Iq/IOOO and p - \. 



small J ^c ojo, the relaxation rates in each channel are equal. 
In contrast, for large AOC J > ojq, the four relaxation rates 
become substantially different. As an example, we note that 
the decay from the excited state |1) to the intermediate state 
|3) is much faster than that of |3) to |4) if the AOC is strong 
and if the coupling constant / has the same sign as p. This 
indicates a possibility to trap the system in the intermediate 
state 1 3), which can either be a symmetric or an antisymmet- 
ric state depending on the sign of the AOC constant /, see 
Eq. jTSh . Interestingly, the rate 7|i>^|3> can be comparable to 
7li>^l2>. However, within a practical AOC strength /, this is 
unfeasible for the small separation r2i = /lo/lOOO shown in 
Fig. El 

We note that related effects due to strong coherent driving 
have already in detail been discussed in a number of previ- 
ous works in atomic systems subject to strong external driv- 
ing fieldsSiSi"— , and recently also in superconducting circuits 
strongly driven by a microwave field"^^'^^. 

It is also interesting to discuss the role of pure dephasing in 
this context. In the dressed state basis, the decoherence from 




FIG. 7. (Color online) The dipole-dipole induced shifts QJ [(i) blue 
line], Q~ [(ii) dashed red line] and Qd [(iii) dash-dotted green line] 
as a function of the qubit separation r2i/Ao at a fixed frequency ojq. 

pure dephasing becomes 

^^p = J] ^j{Y^/2, \ejXej\ - \gj}{gj\}p 
j 

=^{(a^-b^)%,\l}{l\-\^m}p 
+ if{4aVr^,|l)<4| + |4)<l|}p 
+ ^{r^,|2)<3| + |3)<2|}p 

|l><4| + |4)<l|}p. (33) 

We find that independent of the AOC strength, the pure de- 
phasing induces incoherent energy exchange between two in- 
termediate states 1 2) and |3). However, it gives rise to inco- 
herent energy exchange between |1) and |4) only if the AOC 
is large. For small AOC the pure dephasing rates of |1) and 
|4) are equal to those of an isolated qubit. This part becomes 
small if the AOC is large. 

Next, we analyze the effect of the dimensionality of the 
environment. In the atomic model of coupled qubits, typi- 
cally a 3D reservoir is considered^. Then, the dipole-dipole 
shifts (DDS) Qj = - Ul^ -\- Q^ j for small distances r2i de- 
cay as r~^ and for larger distances vanish with oscillations 
around zero and an overall decay proportional to r^/ . In our 
model, the reservoir effectively forms a one-dimensional open 
space"^"^. In ID, the DDS Dj is much smaller for small separa- 
tion r2i < /lo/ 10 than in 3D, see Fig.[7l The DDS oj^^^ -oj^^^ 
between the two intermediate levels |2) and |3) quickly decays 
proportional to r'j^ as the separation increases to r2i > Ao/50. 
This shift is about IOO70 for r2i = /lo/lOOO but decreases to 
less than 70 if ^21 < /lo/50. Furthermore, Qj reduces to a si- 
nusoidal oscillation if r2i > Aq/IO, as demonstrated for two 
qubits with coupling mediated by one-dimensional plasmonic 
waveguides^ ^. Thus, the DDS effects on the evolution of sys- 
tem is small except for trivial coherent shifts in the transition 
frequencies. We note that our results is different from the pre- 
diction in^^ in the case of a distance much smaller than the 
wavelength. Our model calculates the principal value integra- 
tion for the DDS, whereas in^\ an approximation involving 
Green's function was used. In our case, the DDS in ID space 




FIG. 8. (Color online) (a) Pictorial representation of a photon emis- 
sion from qubit Ql in an open 2D space, (b) The corresponding 
situation in ID space. The results are obtained by numerically solv- 
ing Maxwell's equations for an oscillating source with the finite- 
difference time-domain method (FDTD)^^. 



exponentially decays. This decay is similar to the 3D case in 
which the loss from diffraction is also negligible if the emis- 
sion from the donor directly points to the acceptor. 

Also the collective relaxation rates of the qubit in ID are es- 
sentially different from the higher dimensional case. In the 3D 
case, the decay cross rate yu tends to 70 for small separations 
^21 ^ ^- For larger distances, it oscillates and overall decays 
rapidly as the separation r2i increases^^. In contrast, in an ID 
bath, this cross damping rate oscillates as a cos function of the 
separation r2i , but does not have an overall decay with increas- 
ing distance, see Eq. ([25]). This can easily be understood by 
noting that in ID, the coupling Hamiltonian essentially does 
not carry distance information, and photons emitted from one 
atom are effectively guided to the other atom via the ID space. 
In contrast, in 3D, photons have a high probability to miss the 
other atom, as shown in a rather pictorial way in Fig.[8l This 
feature is comparable to atoms embedded in a suitable cavity, 
which also can render the environment predominantly one- 
dimensional. 



B. Ground state entanglement 

Creating and protecting high entanglement in a controllable 
way is of interest, as the decoherence arising from the cou- 
pling of a system to its environment tends to rapidly destroy 
entanglement. This is particularly difficult in solid-state su- 
perconducting circuits. 

In the following, we first demonstrate the possibility to cre- 
ate ground state entanglement based on the modification of 
dissipation due to the strong coupling between the qubits. In 
particular, no external coherent driving is required, and we ne- 
glect the corresponding Hamiltonian ///. Thus we will show 
that in our system, the dissipation can be a useful resource for 
quantum entanglement without any coherent driving field^^. 
Note that the usefulness of dissipation has recently also been 
demonstrated in a carefully engineered local environment—. 
The possibility of ground state entanglement in superconduct- 
ing qubits has been discussed for vanishing distances between 



two qubits2^. However, in this study, neither the operation 
time nor the collective coherent driving induced by by dissi- 
pation was considered. 

As discussed in Sec. lIIIAl in the dressed state basis all up- 
per states of our system decay to the ground state, see Fig.[3l 
We have verified this by numerically solving the master equa- 
tion ([3Qb with Hi = 0, as shown in Fig. [3a). To assess the 
ground state entanglement, we calculate the entanglement of 
formation of the two qubits in the magic basis asS 



where 



E(p) = S(C(p)), 



1 + ^^^^ 



S(C) = h 



h(x) = -X log2 X- (I - x) log2(l - x) , 



(34) 

(35a) 
(35b) 



with concurrence C(p) = max{0,Ai - A2 - A3 - A4} = 
max{0, 2Amax - Tr[R]}. The At are the eigenvalues, in decreas- 
ing order, of the Hermitian matrix R given by 



p = {CTy® Cry)p*(Cry (Ty) . 



(36a) 
(36b) 



For small coupling strengths J ^^ ojq, our system reduces to 
the usual configuration studied, e.g., in^'^^. In essence, the 
ground state tends to a product state, and neither entanglement 
nor coherence occurs, see Fig.[9l 

Next, we focus on the case of large AOC, both without pure 
dephasing and with pure dephasing of F^ = O.Olyo- Surpris- 
ingly, for the AOC |/| larger than the resonance frequency of 
the qubits, the ground state tends to a Bell state 



\^±) = —=(\gug2)±\eue2}) 
V2 



(37) 



with large fidelity, such that the entanglement becomes high. 
The stronger the AOC is, the higher is the fidelity of the Bell 
state, and subsequently the larger the entanglement. For a 
large coupling |/| > 5^o, the fidelity of the Bell states ex- 
ceeds 97% in the steady state if no pure dephasing occurs, 
as shown by the solid lines in Fig. [9ta). The ground state 
is close to |0_) for a positive /, while it becomes |0+) if 
/ < 0. By dynamically tuning the sign of the coupling^—, 
the ground state can thus be switched between the two Bell 
states |0+). Alternatively, the bias flux can be inverted around 
the optimal point because this also inverts the sign of the 
AOC, see Fig. [2l As expected, the population in the product 
states 1^1,^2) and |^i, ^2) is even and about 0.5, see Fig.[9tb). 
Related to the entanglement, also near-maximum coherence 
K^i, ^2lpki, ^2)1 ~ 0.5 can be created via the dissipation. The 
dashed lines in Fig. [9] show the influence of pure dephasing. 
If the pure dephasing is low, then it only slightly reduces the 
fidelity of the achieved Bell state. For example, the entan- 
glement for 2.8 < \J\/ojo < 6.3 decreases to about 0.87 cor- 
responding to a 95% fidelity of the states |0+) if a pure de- 
phasing of F^ = 0.0 170 is taken into account. The coherence 
K,^i, ^2lpki, ^2)1 is also only slightly reduced. 




FIG. 9. (Color online) (a) Fidelities of the Bell states |0+> [thin blue 
line (i)], |0_> [green lines (iii)], and entanglement [red lines (ii)] as 
a function of the AOC /. (b) Fidelities of the product states ki, ^2> 
[blue lines (i)], |gi,g2> [green lines (ii)], and coherence [thick cyan 
lines (iii)] as a function of the AOC /. In both subpanels, the pa- 
rameters are /? = 1 and r - /lo/lOOO. Note that the coherence was 
shifted upwards by 0.5. The solid lines indicate the case without 
pure dephasing, while the dashed lines include a pure dephasing of 

r^ = o.oiro. 



Next to the achievable fidelity, also the operation time is 
a crucial characteristic of any quantum gate. To explore the 
temporal dynamics of our system, we start in an initial prod- 
uct ground state \g\,g2) with AOC set to zero. Then we switch 
on the AOC in a switching time tsw After a certain time, we 
switch off' this coupling again. A complication in the theoreti- 
cal analysis arises from the fact that with the time-dependence 
of /, also our transformation to the dressed state basis be- 
comes time-dependent. The bare basis is transformed to the 
dressed-state basis via an unitary operation defined as 



% = 



^ a ^ {) b\ 

fi a 

a -/3 

-b a 



(38) 




0.05 , 0.1 

^0^ 



FIG. 10. (Color online) Time-evolution of entanglement for time- 
dependent couplings between the qubits. The time-dependent cou- 
pling strength / is shown as thin blue line (i). The amount of entan- 
glement is shown as red line (ii). The fidelity of the Bell state |0_) is 
indicated as the dashed green line (iii), and the total excitation from 
the ground state as dot-dashed cyan line (iv). Subpanel (a) shows the 
case of slow switching of the AOC / with switching time tsw == 7o^- 
(b) shows the corresponding results with fast switching tsw = O-OITq ^ 
(b). The results are obtained for r2i = /lo/lOOO and include a pure 
dephasing of F^ = O.Olyo- 



introduces a coherent contribution 

Hx = -in%dt%^ 



(39) 



to the transformation of Hamiltonian^. We thus obtain the 
master equation 

j 

Results of the numerical solution of Eq. (l4Ql) are shown in 
Fig- [El Initially, the system is |4). The AOC is swept from a 
negligible value yo to the large value 5^o and then back to a 
vanishing coupling. In Fig. [I3a), the switching of the AOC 
occurs rather slowly, with switching time on the order of the 

In this case, the system approxi- 



relaxation time, tsw = 



To" 



Taking into account the time dependence of the coefficients mately remains in the time-dependent ground state throughout 
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the whole evolution, with less than one percent of population 
in the three upper states during the switching. Due to the evo- 
lution of the ground state, the system changes from a product 
state 1^1,^2) to a highly entangled state |0_) with fidelity of 
0.98 and corresponding entanglement E{p) = 0.95, and back 
to the product state again when we turn off' the AOC. 

Next, we assume a faster switching time, see Fig. [TOtb). 
Although the switching time is only tsw = O.Oly^^ still more 
than 98% population remain in the ground state |4) at all times. 
The fidelity and entanglement are similar to the case of slow 
switching. Thus, a fast and robust entangling operation can 
be realized in our system with the use of a tunable AOC. Both 
highly entangled states |0+) can be prepared from the product 
state 1^1,^2) via tuning on a large AOC, without any external 
coherent driving fields. During this operation, the system re- 
mains in its ground state. Only non-adiabatic eff'ects during 
faster switchings of the AOC give rise to a partial excitation 
of the upper states. Thus, this operation is nearly free from 
decoherence. 

Interestingly, in Fig. [I3b) it can be seen that the fidelity 
of |0_) oscillates rapidly after the switching sequence is fin- 
ished. These oscillations at constant AOC arise from the slight 
excitation of upper states due to the fast switching of a large 
AOC. This partial excitation induces coherences which oscil- 
late at the transition frequency. The excitation pe of state \E) 
then gives rise to an oscillation of fidelity of |0_) with am- 
plitude ^/pE. For example, the excitation pe is about 0.0025 
after the AOC is switched off' at 70^ = 0.1 1. Correspondingly, 
the fidelity of |0_) oscillates with amplitude 0.05 = VO.0025 
around the mean value 0.5. 

We finally turn to shortest switching times of order of the 
inverse transition frequency of the qubit. Since the AOC / ex- 
ceeds the transition frequency 2ojo of the doubly excited state 
\E) at its maximum value, it should in principle be capable 
of exciting the qubit system to the excited state to a large de- 
gree. Similarly, due to non-adiabatic eff'ects, large residual 
entanglement after the switching should be possible. This is 
demonstrated in Fig.fTTIa), which shows the entanglement be- 
tween the two qubits in dependence on both the AOC strength 
/ and the switching time tsw The duration of the AOC is 2tsw, 
and To = 2n/ojo is the carrier period. We find that the entan- 
glement of the system can be large after the AOC turn off' if 
the area of the AOC pulse satisfies 




0.03 



FIG. 1 1 . (Color online) (a) Amount of entanglement as a function of 
the switching time tsw and the coupling /. The white dashed lines in- 
dicate pulses with area 3:7r/2, Sn/l, In/ 2, 9:7r/2, 1 l;7r/2, 13:7r/2, 15;r/2, 
respectively, (b) shows the time evolution of the entanglement gen- 
eration for a pulsed AOC /. The time evolution of the coupling / is 
shown as blue line (i). The resulting entanglement is plotted as red 
line (ii). The fidelity of the state |0_> is shown as green line (iii); that 
of the state |0+> as cyan line (iv). Finally, the state population of \G) 
is shown as green dashed line (v), and the population of \E) as cyan 
dashed line (vi). Note that the lines (i), (v) and (vi) are shifted down 
by -1 for better clarity. In (b), the switching time is tsw = 0.00277^^ 
The other parameters are: F^ = O.OI70 and r2i = /lo/lOOO. 



2QRabitsw = nn/2 , 



(41) 



where n e {3,5,7,9,...}, as indicated by the dashed lines in 
Fig-EIa). Here, the eff'ective Rabi frequency is given by 



^Rabi - ^0 ' 



[OJO \2tswOJ0 



-1 



(42) 



The high residual entanglement goes along with a high popu- 
lation of the excited state \E). Efficient population of \E} can 
be expected for 2tsw/To ~ 1 . As shown in Fig. fTTT b). the states 
\E} and |G) are then evenly populated after the switching op- 
eration is finished. Interestingly, the entanglement of the sys- 
tem can be E(p) = 0.98 after an operation with / ~ l.Scoo and 
2tsw ~ 7^0/4. After the switching, the population oscillates 



between the two Bell states |0+), as they are not the eigen- 
states of the system. The entanglement itself, however, does 
not show this oscillation. 

We thus have demonstrated an experimentally realizable 
mechanism for the generation of entanglement and for gen- 
erating maximum coherence in time scales of order of the in- 
verse qubit transition frequency, exclusively driven by dissipa- 
tion. In contrast to intuition, the dissipation turns into a favor- 
able resource if combined with the strong coupling between 
the two qubits. It can be seen from Fig.[9]that the ground state 
is continuously tunable with the restriction that the population 
of the excited state \E) remains smaller than 50%. 
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FIG. 12. (Color online) Level diagram and relevant coherent and 
incoherent processes for the generation of stationary state entangle- 
ment with a single continuous wave driving field. As an example, 
we use the parameters correspond to those indicated by the point P 
in Fig. [Ha). 




FIG. 13. (Color online) Optimizing the entanglement generation 
with a single cw field. The figure shows the fidelity of the Bell state 
|^_) as a function of the detuning 6. The solid lines are drawn for 
a distance r2i = /lo/50, the dashed lines show a smaller separation 
^21 = /lo/lOOO. The two upper blue lines correspond to the ideal sys- 
tem without pure dephasing, whereas two lower red lines shows the 
results taking into account a pure dephasing F^ = O.Olyo- 



C. Excited state entanglement 

In the previous section, we have considered entanglement 
generation in the ground state. Here, we analyze the genera- 
tion of entanglement between two excited qubit states. Usu- 
ally, such excited states are subject to dissipation, and thus 
have only a limited life time. Therefore, in the spirit of deco- 
herence free subspaces2ii22, we will analyze the possibility of 
generating trapped entangled excited states which are immune 
to the dissipation of system, and focus on the states 




Y2 



(43) 



FIG. 14. (Color online) (a) Steady-state entanglement with a sin- 
gle continuous-wave driving field as a function of the driving field 
Rabi frequency Q and the coupling J. The upper panel (a) shows 
the results for parallel dipole moments, p = 1, of the two qubits. 
The lower panel (b) shows the corresponding results for anti-parallel 
dipole moments, p = -1. In both cases, the frequency of the driv- 
ing field is chosen as half the transition transition frequency of the 
|1) -^ |4) transition. Thus, for large coupling /, all transitions are far 
off'-resonant with the field. The phase is chosen as 0i2 = 0. 



In atomic systems with weak coupling, it is challenging to 
highly populate the states in Eq. (|43l) . The anti- symmetric 
state is decoupled from the driving field dynamics if both 
atoms are subject to the same driving field with the same 
phase, as it is usually the case for two atoms with a distance 
of order of the transition wavelength^^ . On the other hand, it 
is also difl&cult to robustly create the symmetric state because 
it couples to both the collective ground and excited states^. 

Superconducting circuits diflfer from atomic systems in 
many ways, and these diflTerences can be exploited to popu- 
late the states Eq. (|43]) . Methods to robustly prepare the an- 
tisymmetric state have been discussed in our previous work^^- 
using adiabatic passage. For this, first a diflTerence between 
the two transition frequencies of the qubits is introduced dur- 
ing the population stage, and then this diflference is tuned to 
zero to protect the prepared state. Accordingly, the generated 
antisymmetric state decays at a very small rate. Ojanen et. al^ 
have also suggested a superconducting circuit to prepare the 



M9^ 




FIG. 15. (Color online) Interpretation of the steady state entangle- 
ment shown in Fig. O in terms of the population of the magic state 
basis states. The four sub-panels show the fidelity of the states (a) 
|0_), (b) |^_>, (c) |0+>, and (d) |^+>. The figure is drawn for parallel 
dipole moments (p = 1). 



antisymmetric state. They designed the circuit geometrically 
such that a microwave n pulse inserted in a coplanar super- 
conducting cavity anti- symmetrically drives the two coupled 
qubits to their antisymmetric state. However, in this setup, 
also the environmental noise in the bath (cavity) couples to 
the antisymmetric state. As a result, the antisymmetric state 
will decay fast. Furthermore, the geometry is fixed and the 
transitions involving the symmetric state are forbidden. Thus 
it is difl&cult to highly populate the symmetric state. 

Motivated by this, in the following, we propose two meth- 
ods to prepare and trap the entangled anti- symmetric and sym- 
metric states. 



D. Trapping assisted by a single field 

Here, we discuss the creation of stationary entanglement by 
using only one continuous-wave driving field. As both states 
|^+) are maximally entangled, the system acquires large sta- 
tionary entanglement if either of these two states is trapped 
with high fidelity. 

The considered configuration is shown in Fig. [121 The sys- 
tem is driven by one cw field with Rabi frequency Q.. In the 
interaction picture defined by the unitary transformation 





(44a) 


ojRi = 2(jJi -\- (jj\ , 


(44b) 


<^R1 = <^L + <^4 , 


(44c) 


^/?3 = ^L + ^4 , 


(44d) 


(j)M - OJa , 


(44e) 



the master equation for the system density matrix p takes the 
form 

P = ^[^0 + ^/ + Hbie,p\ + ^spp + ^TPP , (45) 
n 
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where ^o = Yjj^j\j){j\ is the free-evolution Hamiltonian. 
The interaction Hamiltonian is 

Hi = no. [a [a + jSpe-'"^'' ) + ^ (y^ + ape-'"^'' )} Ru 
+ no. [a [a - Ppe-'"^^' )-^b{p- ape-'"^^' )} 7?34 
+ no. [a {p + ape-'"^'' ) - /? (a + jSpe-'"^'' )} R24 
- no. [a {p - ape-'"^'' )-b{a- fipe-'"^'' )} 7?i3 -h H.c. , 



and the detuning are defined as 

Ai =26, 

A2 = S -\- <jJ2- 

A3 = S -\- (JJ3 - 
A4 =0, 



OJl + <jJ4 

2 

OJl + OJ4 



(46) 

(47a) 
(47b) 

(47c) 
(47d) 



where S = (€01-0)4) /2-ojl. The detuning S can be controlled 
via the frequency of the applied TDMF. The phase 02 1 = ^l^21 
of the driving field depends on the frequency of the driving 
field and is non-zero. Therefore, the channel relevant to |3) 
needs to be included. Note that the coherent part Hbie result- 
ing from DDI remains unchanged under the applied transfor- 
mation. 

To analyze the system performance, we choose a driving 
field with Rabi frequency Q = 8O70, and an AOC |/| = 5^o- 
Sweeping the frequency of the driving field, we find an opti- 
mal detuning Sopt = 2.5yo leading to maximum fidelity in the 
target state, see Fig. [131 In the absence of pure dephasing, one 
can create the Bell state |^_) with a high fidelity in a broad fre- 
quency window for both considered separations r2i = Ao/50 
or r2i = /lo/lOOO. At optimal detuning S = 2.5yo, the fidelity 
is one. Including pure dephasing with F^ = 0.01 70. the fidelity 
reduces to a small value of about 0.22 for small separation. In 
contrast, one can still trap the state |^_) with a fidelity of 0.83 
for larger distances r2i = Ao/50. The reason is that in the case 
of small separation, the decay from |1) to |3) and thus to |^-) 
is several orders of magnitude smaller than that to |2), which 
decays quickly to the ground state |4). But for r2i = /lo/50, the 
decay rate 7|i>^|3> is comparable to 7|i>^|2> and several orders 
of magnitude larger than 713)^14). Therefore, most population 
decays to state |3) for larger separation. Similarly, one can ob- 
tain the comparable trapping in the state |^+) at ^ = 2.5^0 if 
/ = -SojQ. 

Figure[T4l shows the entanglement of the system against the 
driving field Rabi frequency and the AOC strength in the ab- 
sence of pure dephasing. The system is highly entangled in 
three regions (I, II, III) in either case of parallel (p = 1) or 
anti-parallel (p = -1) dipole moments. In these regions the 
AOC strength / is larger than the transition frequencies. In 
region IV, the system is weakly entangled, and in regions V, 
VI and VII, it is almost disentangled. 

The structure of entanglement shown in Fig. [I^a) can be 
understood as follows. We decompose the system state into 
the the magic state basis^, show the fidelity of the four magic 
states in Fig. \T5\ and discuss the diflTerent regions separately 
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FIG. 16. (Color online) Steady-state entanglement and interpreta- 
tion in terms of Bell state populations against the driving field Rabi 
frequency Q and the coupling /. The figure is similar to Fig. [141 but 
includes pure dephasing F^ = 0.0 lyo- (a) shows the amount of entan- 
glement, (b) the fidelity of the Bell state |0_>, and (c) that of the Bell 
state |^_>. The results are shown for anti-parallel dipole moments 
(P = -1). 



in the following. In region I, the driving field is weak, and oflT- 
resonant from all transitions. As a result, almost all population 
is in the ground state, which is close to the Bell state |0_) for 
the considered large AOC /. This expectation is confirmed in 
Fig.[T5ta). From Fig. \T5\h) we see that interestingly, the an- 
tisymmetric state can be prepared using a single driving field 



in region II. The field first pumps the system into the interme- 
diate state |2) and then into the excited state |1), which subse- 
quently decays to the ground state via three channels: the sym- 
metric, the antisymmetric and the two-photon channel. Essen- 
tially diflferent from the weak AOC case, the decay from the 
antisymmetric state to the ground state is strongly modified to 
be much smaller than other channels by the strong AOC, as 
shown in Fig. \5\ Thus eflTectively the system is pumped into 
the antisymmetric state and remains trapped in it. This leads 
to high entanglement. For example, E > 0.95 is achieved for 
Qi = Q2 = 8O70 and / = Sojq, see Fig. [Hta). It should be 
noted that this state is created in the steady state of system. 
Therefore it is independent of the initial state protected from 
the relaxation, thus forming a decoherence-free subspace. For 
negative /, in region III, almost all population decays to the 
ground state which is now close to the Bell state |0+), see 
[T5lc). The reason is that the decay from |1) to the intermedi- 
ate states |3) and those of these intermediate states are equal 
but much larger than that to |2), see Fig. [5] Thus only a small 
amount of population resides in the symmetric state |3) which 
means that the fidelity of the symmetric state is small. Again, 
the system is highly entangled in steady state. At certain driv- 
ing fields, e.g., in the region IV, the upper states |2) and |3) 
are partly excited. The system then dominantly is in a super- 
position of |0+) and |^+) with equal weight. This results in 
an unentangled state. Similarly, the system is unentangled in 
region VI. In region V, the AOC is small. This reproduces the 
regular case of two weakly coupled two-level systems. The 
system is weakly excited by the oflT-resonant driving field, and 
the population of all four states are similar, see \T5\ Thus the 
system is unentangled. 

Further investigation shows similar entanglement proper- 
ties as a function of the AOC and driving field for opposite 
dipole moments Ji ][ J2, see Fig. [T4Tb). In this case, the 
symmetric state |^+) can be trapped by a strong driving field 
in the region of negative coupling. 

If a pure dephasing of F^ = 0.0 lyo is included, again the 
entanglement properties shown in Fig. [T6ta) depend on the 
eflfective distance of the qubits. For small separation such 
as r2i = /lo/lOOO, the system is entangled only if almost 
all population is in its ground state, see Fig. [16] (b). If the 
AOC is small, populations in all three upper states are nearly 
equal and have similar structures as the state |^_) shown 
in Fig. [T6tc). Their maxima all are about 0.26. As a re- 
sult, the entanglement is nearly zero in the regions marked 
as y, VI, VII in Fig. [T6ta). For small separations, either the 
decay rate 7|i>^|3> or 7|i>^|2> is orders of magnitude smaller 
than 7o, see Fig. (6] The system thus reduces to a ladder-type 
three level system. In addition, the incoherent exchange due 
to large pure dephasing reduces the diflTerence in population 
of states |2) and |3). Thus, the entanglement and all fidelities 
of the eigenstates becomes symmetric with respect to / = 0. 
The fidelity of state |^_) is shown in Fig. [T6tc), and the fi- 
delity of states |1) and |2) have a similar structure. In regions 
V, VI and VII, all eigenstates are evenly populated, such that 
the entanglement vanishes. 

For larger separations such as r2i = /lo/50, an entanglement 
of 0.58 can be observed in region II of Fig. [TTta). It arises 
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FIG. 17. (Color online) Steady-state entanglement and interpreta- 
tion in terms of Bell state populations against the driving field Rabi 
frequency Q and the coupling /. The figure is similar to Fig. [161 
but generated for larger qubit separation r2i = Ao/50. (a) shows the 
amount of entanglement, and (b) the fidelity of the state |^_). 



mainly from contributions of the state |^-) because it is still 
created with a high fidelity of about 0.83, see Fig. fTTl b). 

We conclude by analyzing the time evolution of the system 
driven by one field for large separation r2i = /lo/50, as shown 
in Fig.[T8l Interestingly, the population first is pumped to the 
upper states |2) and |1), which decay to |3). From time yot = 3 
on, the damping of 1 1 ) is dominant. Therefore, the population 
of state |3) continually increases to is maximum value of 0.83 
at about yot = 80. The system then remains in this steady state 
with E(p) = 0.58. 



E. Trapping assisted by dark-state techniques 

Coherent population trapping (CPT) and the STIRAP tech- 
nique are a convenient tool to prepare atoms in superposition 
states without decoherence throughout the preparation. Typ- 
ically, two driving fields are applied to an atom to generate 
a A-type level scheme. If the driving fields are chosen suit- 
ably, the system is driven into a dark state, which is a coherent 
superposition of the two ground states, with relative weights 
depending on the driving field parameters. Here, in contrast 
to previous implementations, we discuss coherent population 




200 



FIG. 18. (Color online) Time evolution of steady-state entangle- 
ment generation with a single cw field. The red curve shows (ii) 
the amount of entanglement, the blue thin line (i) shows the fidelity 
of state |^_>, the green dashed line (iii) that of |0_>, the cyan line 
(iv) that of |^+), and the orange dot-dashed line (v) that of |0+>. The 
parameters are r2i = /to/50, /? = 1, (5 = 2.57o. T^ = O.Olyo and 
Q = SOyo- 




FIG. 19. (Color online) Level diagram and relevant coherent and 
incoherent processes for entanglement generation using dark state 
techniques. 



trapping and adiabatic population transfer between two entan- 
gled states. Thereby, the entangled states are protected from 
decoherence, and the preparation is robust against noise in the 
field. 

To enable dark state techniques, we apply two driving fields 
denoted as control field Dc and pumping field flp. These fields 
drive the two-photon and the |1) ^ |3) channels, respectively, 
as shown in Fig [191 The transitions |3) ^ |4) and |1) ^ |2) 
are far oflT-resonant and can be neglected. The control field has 
frequency oj^^-cdr^-ojc, the two-photon channel is driven 
ojp. We define a unitary transformation 



with frequency co]^ 



"-L' 



JCORjt)\j}{j\ 



(jjRi =0J4-\-ajp, 



(48a) 
(48b) 
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(x)R2 =(x)A+(x)c, 


(48c) 


iOR3 = (j}^ + (j}p - 5c , 


(48d) 


com = C04 ■ 


(48e) 



The one-photon detuning is defined as ^^ = (oji - 0)4) - ojp 
and 6c = (oji - 0J3) - ojc. This gives 



Ai=Sp, 




(49a) 


A2 = ((JJ2 - CO4) - (COI 


-(jJ3)-\-Sc, 


(49b) 


A3 =Sp-Sc, 




(49c) 


A4 = 0. 




(49d) 



In the RWA, the master equation for this scheme reads 



Z^I^'X^'I' 



-[^/,p] + ^p + ^rpp, (50) 
n 



where 



Hi = n^c [a (p + ape-''^'^) - h (a + Ppe-''^'^)] R24 
- n^c [a [p - ape-'"^^) -b{a - Pp''"^')] R13 
+ 2nap(a^ - b^)[l + e-''^p)Ri4 + H.c. . (51) 



Here, we have neglected the off-resonant excitations |1) ^ |2) 
and |3) <-^ |4), as well as rapidly oscillating contributions from 
Hbie- We also exclude rapidly oscillating parts relevant to 
y?2\^^ and y?^)^^ in -^P yielding ^p. The phases are cpj = 
27i{r2\ I Xq){(jl> j I cjJq) with j e {c,p}. Unavoidable, the state |2) 
is coupled to the ground state by Qc- 

In Fig. [20I (a), we demonstrate a coherent trapping of the 
intermediate state |3) via the dark state approach. This state 
is of interest since it can approach the maximally entangled 
Bell states |^+), and the decay rate can be modified to be very 
small by the strong AOC. The optimal trapping is obtained for 
Sp = Sc = 0, flp = 4yo, and Q^ = To- Then, the fidelity of 
|^_) is about 0.87 corresponding to an entanglement of 0.66. 
Only about three percent population remains in the state |1). 
The intermediate state |2) is populated to about 6%. Note that 
the state |3) could be either the symmetric state for / > or 
the antisymmetric state for a negative /. The population in 
the state |3) can be controlled by the field flc, as shown in 
Fig.[2Qtb). For small Q^, the state |3) is highly trapped, as the 
corresponding dark state has a high contribution of |3). For 
Qc » ^/7, correspondingly the system is trapped in |4) and 
1 2). When the pumping field flp becomes comparable to the 
controlling field Q^, the two entangled states are equally popu- 
lated. The system is in a superposition state of two maximally 
entangled states, which leads to vanishing overall entangle- 
ment as the superposition state is a separable state. As the 
controlling field increases further, the state |2) is excited. The 
population in this state is similar to that of the ground state. 
These two states together form a separable mixed state, and 
the system becomes disentangled. 

We also investigated the time evolution in the dark- state ap- 
proach. Similar to the trapping assisted by a single field, we 
found that it takes a rather long time (> 60^^^) to reach the 
maximal value of entanglement and fidelity of the target state. 




p 



FIG. 20. (Color online) Entanglement generation with the dark state 
technique. In (a), the results are shown against the detuning 6 p. The 
red curve (ii) shows the amount of entanglement, the blue thin curve 
(i) the fidelity of |^_>, the green dashed line (iii) that of |4>, the cyan 
line (iv) that if |2>, and the orange dot-dashed line (v) that of |1>. The 
parameters are Sc = and Qp = 70 • Corresponding results shown 
against the control field Rabi frequency Q^- The parameters are Qp = 
470 and 6p = Sc = 0. In both subpanels, the other parameters are 
rii - /^o/50, r^ = O.OI70, 7o = 0.0270, / = 5coq and p = \. 



The STIRAP technique^^ is another robust method to co- 
herently prepare a target state. In the following, we apply 
STIRAP to our system. For this, we use the same configura- 
tion shown in Fig.[20land sech-shaped pulses defined as 

Qc = 3flosech[7i(t - 5tc)/2tc] , (52a) 

Qp = Qosech[7r(^ - Szp - Td)/2Tp] , (52b) 

where t/ with / g {c, p} are the durations of the pulses, and 
Td the delay between the pulses. The factor 3 is introduced to 
balance the prefactors of 7?i3 and Ru in Eq. (ISTI) . 

We found that the CPT is optimal in two-photon resonance 
Sp = 6c = for Tp = Tc = T. As shown in Fig. [211 using 
parameters yoT = 5, tj = 3t and Qq = IO70, about 70% pop- 
ulation is transferred from the ground state |0_) to the state 
|3), which tends to |^_) for positive / or |^+) if / < 0. 
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FIG. 21. (Color online) Coherent population transfer from the 
ground state to state |3> via the STIRAP technique. The Stokes pulse 
Qc (blue thin dashed line) precedes the pumping pulse Qp (red thin 
dashed line) by ISy^^. In the figure, both pulses are shown nor- 
malized. The other curves show the amount of entanglement [red 
line (ii)], and the fidelity of states |^_> [blue thin line (i)], |4> [green 
dashed line (iii)], |2> [cyan line (iv)], and |1> [orange dot-dashed line 
(v)]. The durations of both pulses are chosen as r = ^Tq^- The detun- 



ings are Sp 



■■ 6. The other parameters are / = Scdq, r2\ - /lo/50, 



r^ = O.Olyo, To ^ 0.0270 and /? = 1. 



IV. SUMMARY 

In this work, we have derived and discussed the master 
equation for two interacting qubits coupled to a common ID 
bath. In particular, we have focused on the case of strong 
interaction between the qubits, with coupling rate exceeding 
the respective transition frequencies of the two qubits. Deco- 
herence of the qubits is taken into account by coupling the 
system to a bath of harmonic oscillators. We have found 
that the strong coupling between the two qubits can lead to 
ground state entanglement even in the absence of any driving 
field, with the entangled state populated by spontaneous emis- 
sion. We have also analyzed the case of a time-dependent 
strong qubit-qubit coupling. We found that starting from an 
initially unentangled ground state of the two qubits, a dy- 
namical switching of the qubit coupling can lead to entan- 
glement with high fidelity on time scales on the order of the 
inverse qubit transition frequency. Employing a single exter- 
nal driving field, both symmetric and antisymmetric collective 
qubit states can be prepared and maintained with high fidelity, 
thus also leading to entanglement. Finally, using two external 
fields, we also demonstrated the possibility for entangled state 
preparation via a dark state approach or STIRAP technique 
in our configuration. Possible applications of our scheme in- 
clude nonlocal quantum gates and quantum memories. 
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Appendix A: Decay rates 

In this appendix we derive the cross decay rate and the de- 
cay rate of isolated qubits. We model the bath as an open 
ID transmission line-i with total inductance Ly = LI and ca- 
pacitance Cr = Lc, where / and c are the inductance and ca- 
pacitance per unit length, respectively. L is the length of the 
line. The zero point fluctuations of the current in the trans- 
mission line is Ir = ^/fiaJr/2LL The coupling 77^^ introduced 
in Eq. ([20l) reads^ 



H 



,(0- 



ek* 



^'^Mpflr. 



(Al) 



The single photon interaction Hamiltonian between the 
qubits and the bath takes the form 



Hbsp = ^ X ^^'^ Z "^^k [^'''^k + H.c) , (A2) 

/=i k 

with polarization ek of the current, which only takes two direc- 
tions, r^^^ indicates the position of the /th qubit in the vacuum 
field. The wave number of the ^th mode is /: = ^r/v, with the 
wave phase velocity v = 1 / ^[lc. 

According to Eq. ([2Ql) , the cross decay rate is given by 

721 = ^ 7id{oj- ^)\rikf' (ek • ^ cos(^r2i) , (A3) 



where r2i = ^2 - ri. We consider the transmission line as 
an Ohmic bath with a high cutoflf frequency cocut such that 

Yjki^k* ^ ^ 2^^ f!'^ e~^^^'"'doj. Here we have introduced 
an exponential cutoflf^^. Then, the rate yu can be evaluated as 



721 = — sm^ — -^ 



7tS(co- g)— cos(kr2i)e-'^^'^''"daj 

2Lr 



= ro(e)cos(-r2ij , 



(A4) 



where e is the transition frequency between the two involved 
states. For simplicity, we have assumed di = di = a. Here we 
define the decay rate for the isolated qubit as 
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2 {MbIp) 

2Hz 



(A5) 



Similarly, we derived the two-photon decay rates fo and y^. 
Following the diagonalization, the rates in Eq. (l22b are given 
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by 



thermore, 



rii>^l2> = 2(1 + Aabap) [y^-^ + 2a/3py^;^) , 
ri2)^|4) = 2(1 - 4aba/3) (y*,^' + 2ai3p/^>) , 

(cross) 
y\2} 



(A6a) 
(A6b) 



(a' - b') [2aP(y^;^ + y^^^) + p(y[-2^ + r^?)] 



rii>^l3> = 2(1 - 4aZpa;S) {y^;^ - 2a/3py[^^) , 
ri3>^|4> = 2(1 + 4aZpa;S) (r^-^ - 2a/3py[-^) , 

^|3> 



(A6c) 
(A6d) 
(A6e) 



(A6f) 



where the dipole-dipole related contributions are listed in Ap- 
pendix |B] 



(±) w±\J\ 
To = ^ — ro(^o) , 



^0 

.(±) _ ^,(±) . 



_ _ w±\J\ 
ri2 = To COS I 27ir2i/Ao 

2w 

ro = — ro(^o) , 

/2w \ 

ri2 = ro cos — 27rr2i//lo • 
\^o / 



(B3a) 
(B3b) 
(B3c) 
(B3d) 



Appendix C: Dipole-Dipole shifts and bath-induced excitations 



Appendix B: Dipole-dipole induced couplings 

The coupHng coefficients related to the dipole-dipole inter- 
action can be evaluated to give 

^(d,±) = - (^''^,±) + ^(~^,±)) 



,(±) 



71 L 

2 1 



-171 e- + 



^+ ^ ro(w+|/|) ^ p CD 

^^'-^ 71{W±\J\) Jo 0J + {W± 



\J\) 



X COS 



e "^cut doj 



y^^^ r 1 

= -^^9ie ^-^-r(0, -s±) + — 
71 Y S + 

^_ ^ yo(w±\J\) ^ p gj 

^^'^> ;r(w±|/|) Jo oj-(w±\J\) 

X cos I — r2i \e'^cut dco 



where 



(Bla) 



(Bib) 





TT 


/* 


no, S±) + 7:ie' 


+ 


1 


(Blc) 


Q« 


TT 




' 


(Bid) 


^i-' 


71 


e^Y(0,S)-^7iie^ - -r 


' 


(Ble) 




(XJQ 


1 . ^-21 

-F/2;r— 


, 


(B2a) 




s± = - 


w ± 

0J{ 


3 


-h/27r 

^c 


rii' 


' 


(B2b) 



with gc = ^c/^o, and coc is the exponential cutoff frequency 
of the Ohmic bath. T is the incomplete gamma function. Fur- 



The coherent part Hdd includes the dipole-dipole shifts 
Hdds = ^ Tjj ^T^^^\j){j\ in the transition frequencies, where 



(DDS) 

^1 



(DDS) 



-2(a/3 -h ab)pQr^^_^ -h 2{ap - ab)pQr^^ 



+) 



- 2{a^ - b^fOr. , 



(DDS) 



-2(a/3 + ab)p^l^_^ - 2{ap - ab)pQr^^^^ , 
«^3 ~ 2(^/5 - ab)pD.'l^j_^^ -\- 2(al3 -\- ab)pD.'^^_^ , 

JDDS) ^ _2^^^ _ ab)pal,^^ + 2(ayS + aZ7)/7Q^^ _) 
- 2(a^ - b^fa-"^ . 

The bath-induced excitation Hbie is given by 



riBiE — J^DTi7 + J^ 



^BIE 



BIB •■ 



(Cla) 
(Clb) 
(Clc) 

(Cld) 



(C2) 



where 



H 



(1) 

5/£ 



//' 



(2) 



- /(a^ - b^)p (Q^^^) - Q(^^ _^) (7?24P^21 - ^12P^42) , 

- Ka" - b^)p (q-^^) - Q-^ _)) {R^2pRi2 - RiipRiA) , 

- i{a^ - b^)p {^l^_^ - Q^^^)) (R34PR31 - R13PR43) , 

- i(a^ - b^)p (Q-^ _^ - Q-^^)) (7?43P^i3 - 7?3ip7?34) . 



Hbie is negligible because it arises from small dipole-dipole 
shifts and is suppressed by a factor of a^ - b^. 



Appendix D: Decay and pumping in the collective basis 

Here we list the decay, pumping and coherent interaction 
contributions in the collective state basis. They are derived 
from the dissipation in the dressed-state basis. 
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The relaxation in terms of the collective state basis becomes 



&ip = a y\i)^\2) {RiepRei - :zReeP ■ 



\pRee] 



2 11 

+ ^ 7|2>^|4> i RcipRiG - ^jRiip - ijpRii 

+ a^y%'''^ [RiEpRiG + RgipRei] , (Dla) 

^ 2 { 1 1 

^2p = CI y\i)^\2>) {R?>ePRe?> - :zReeP - :zPRee 

2 \ 11 

+ ^ 7|3>^|4> { RodpRdo - :zR33P - :zpR33 

- a^Ti^;'"'^ [Rg3PRe3 + R3EPR3G} , (Dlb) 

^ z . 11 

^TPp = a yTp\ RqePReg - -jReeP - :zPRee 

+ abyrp \ RggPRqg - ^RgeP - -pReg \ - (Die) 



The pumping contributions are given by 



b y\i)^\2) {R2gpRg2 - -RggP - -j^pRgg 

? ( 11 

+ b y\2)^\4) <Re2PR2e - -R22P - :zpR22 

- b^y%'''^ [Re2pRg2 + R2GPR2E} , 
^2P = b^y\i)^\3) \R3gPRg3 - :^RggP - :^pRgg 



(D2a) 



+ b r|3>^|4> i Re3PR3E - -jR33p " tP^33 
+ Z^VLT"^ {^£3P^G3 + ^3GP^3£} , 



^\3) 



(D2b) 



g^TPp = b yTP< RegpRge - -RggP - izPRgg 



+ abyTp \ ReepRee - :zRegP - :^pRge \ • (D2c) 

The dissipation in the dressed state representation introduces 
a coherent contribution in the collective basis as 



^ip =aby\i)^\2) \R2ePRg2 - :zRgeP - :zPRge 



+ aby\i)^\2) \R2gpRe2 - -RegP - :zPReg 



- aby\2)^\4) {Rg2PR2e + Re2PR2g} 

(crc 

M2> 
(crc 

M2> 



- aZpy^^T^^^ [Re2PRe2 + R2EPR2E} 



+ aby\;:'''^ [Rg2PRg2 + R2GPR2G} 



(D3a) 



^2P =aby\i)^\2,) \R3epRg3 - :zRgeP - :zPRge 



I 1 1 

+ aby\i)^\'i) <R3gpRe3 - :zRegP - :zPReg 

- aZ77|3)^|4> {Re3PR3G + Rg3PR3e} 

- aby^^'""''^ {Rg3PRg3 + R3GPR3G} 

+ aby^^'""''^ [Re3PRe3 + R3EPR3E} • 



(D3b) 



Note that the relation a^ > b^ holds. Thus, the decay rates are 
always larger than the corresponding pumping rates. 
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